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Observation of half-quantum vortices in topological superfluid 3 He 
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One of the most sought-after objects in topological quantum matter systems is a vortex carrying 
half a quantum of circulation. They were originally predicted to exist in superfluid 3 He-A, but 
have never been resolved there. Here we report an observation of half-quantum vortices (HQVs) in 
the polar phase of superfluid 3 He. The vortices are created with rotation or by the Kibble-Zurek 
mechanism and identified based on their nuclear magnetic resonance signature. This discovery 
provides a pathway for studies of unpaired Majorana modes bound to the HQV cores in the polar- 
distorted A phase. 


Topology reveals unified features in quantum-matter 
systems as diverse as the chiral A-phase of superfluid 3 He 
at temperatures T ~ 10 -3 K and the vacuum of the Stan¬ 
dard Model at T ~ 10 15 K [1], Those common topologi¬ 
cal properties are manifested in supported topological de¬ 
fects (like vortices, skyrmions etc.) and gapless quasipar¬ 
ticle states. Superfluid phases of 3 He provide a versatile 
platform for studying topological properties of quantum 
matter [1, 2]. One of the most striking consequences of 
non-trivial topology is a vortex which carries half a quan¬ 
tum of circulation in a superfluid or of magnetic flux in a 
superconductor. Such vortices are also predicted to host 
core-bound unpaired Majorana modes [3-5]. The search 
for half-quantum vortices (HQVs) began in chiral Weyl 
superfluid 3 He-A in the late 70’s. Despite promising the¬ 
oretical predictions, HQVs in 3 He-A remained elusive in 
experiments [6-10], as the spin-orbit interaction makes 
them unstable [11, 12]. In the meantime observations 
of half-quantum vortices and fluxes have been reported 
on the grain boundaries of d -wave cuprate superconduc¬ 
tors [13], in chiral superconductor rings [14], and in Bose 
condensates [15, 16]. In Bose systems, however, vortex- 
core-bound fermionic states do not exist, while in super¬ 
conductors only half-quantum fluxes not associated with 
a core (in Josephson vortices or rings) have so far been 
studied. Here we report the discovery of half-quantum 
vortices (HQVs) in the recently found polar phase of su¬ 
perfluid 3 He [17, 18]. 

It is established that 3 He becomes superfluid by 
Cooper pairing in a state with spin 1 and orbital mo¬ 
mentum 1, which is described by a 3 x 3 order parameter 
matrix d M -, where indices p, and j represent the spin 
and orbital degrees of freedom. The polar phase appears 
below the critical temperature T c ~ 10 -3 K when liquid 
3 He is confined in nafen [18], a commercially produced 
solid nanomaterial consisting of long strands aligned in 
the same direction z (Fig. la). In the polar phase the 


order parameter has the form 

A^j = A e^d^rhj. (1) 

Here A is the maximum gap in the quasiparticle energy 
spectrum, <j> is the superfluid phase, d is a unit vector 
of spin anisotropy, and m is that of orbital anisotropy. 
Magnetic field H > 3mT fixes d = icosa(r) + jsina(r), 
where i and j are mutually orthogonal unit vectors in the 
plane normal to H. In the nafen m is pinned parallel to 
the strands, rh || z. The spin-orbit interaction F so oc (d • 
rh ) 2 then affects orientation of d in such a way that the 
distribution of a is governed by the Sine-Gordon equation 
[19] 

V 2 a = (sin 2 fj. / 2£f>) sin 2a, (2) 

where ~ 10 /im is the dipole length and /i is the angle 
of the magnetic field with respect to z (Fig. la). 

A quantized vortex is a linear topological defect in the 
order parameter field which traps non-zero winding <t> —> 
0 + 27rz/ of the order parameter phase 0 on a path around 
the vortex core. Winding of the d orientation, given by 
the angle a, is also possible. Thus, three different types 
of vortices may exist (Fig. lb): single-quantum vortices 
(SQV), spin-current vortices (SCV), and HQVs. SQVs 
are usual phase vortices with v = 1, associated with one 
quantum k = h/ 2 m 3 of mass circulation with velocity 
v s = (ft/27713)V<(> (m3 is 3 He atom mass). In a SCV a —> 
a + 27t on a path around the core. It carries quantized 
circulation of spin current, created by reorientation of the 
d field. For a HQV v = 1/2, and it carries only half a 
quantum k/2 of mass circulation. The apparent phase 
jump by 7r is compensated by the change a —> a + -jt, 
which makes the order parameter single valued. 

Reorientation of d outside of SCV and HQV cores is 
governed by well-known solitonic solutions of Eq. (2). 
The SCV becomes a termination line of two and HQV of 
one 7r soliton in the d field. The soliton connects a HQV 
to another HQV with opposite d winding (SM, Fig. SI), 
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FIG. 1. (a) Sketch of the sample: Cubic container made from epoxy (light blue) is filled with the nafen (blue cube). The 
nafen strands are oriented along the vertical z direction, their diameter is d\ ss 9nm, and their average separation d 2 « 35 nm. 
The space between strands is filled with liquid 3 He. The sample is surrounded by NMR coils made of copper wire (brown 
rectangles). The magnetic field H is applied transverse to the NMR coil axis at an arbitrary angle p from the direction of 
the orbital anisotropy vector m, which is pinned along the nafen strands. The sample can be rotated around the vertical axis 
with the angular velocity Q up to 3 rad/s. (b) Vortex types in the polar phase: The order parameter phase (f> is shown by 
the background colour. The spin vector d is locked to the plane (green disks) perpendicular to the magnetic field H. Within 
this plane hatd rotates by 7r around the HQV core and by 2n around the SCV core. In a tilted magnetic field (p > 0) this 
winding is concentrated in one d soliton terminating at the HQV core, or in two solitons terminating at the SCV core. The 
approximate soliton extents are illustrated by dashed lines. The nafen strands, the m field, and the vortex lines are orthogonal 
to the plane of the picture. The SQV and HQV have hard cores (red discs) of the size of coherence length ~ 40 nm, while the 
SCV has only a soft core (blue disc) of much larger size ~ 10 pm. 


to a SCV or terminates at the sample boundary. The 
soliton width is ~ ^u/sin/r. In the absence of magnetic 
field or when the field is oriented along the nafen strands 
(p, = 0), the spin-orbit energy F so is invariantly at the 
minimum for any d _L m, and solitons are not created. 

In our experiments we use continuous-wave nuclear 
magnetic resonance (NMR). In the superfluid state the 
spin-orbit coupling provides a torque acting on the pro¬ 
cessing magnetisation, which leads to a shift of the pre¬ 
cession frequency from the Larmor value ui l = | 7|IT, 
where 7 = -2.04 • 10 8 s^T" 1 is the gyromagnetic ratio 
of 3 He. The sample regions where spin-orbit energy is at 
a minimum form the main peak in the NMR spectrum 
at the frequency w ma i n with the shift [20] 
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Here flp is the pressure- and temperature-dependent 
Leggett frequency in the polar phase, which characterises 
the spin-orbit torque. 

Within soft cores of topological objects, such as the d 
solitons, the spin-orbit energy is not at minimum, which 
provides a trapping potential for standing spin waves. 
Excitation of these waves leads to a satellite peak in the 
NMR spectrum at the frequency w sa t with the shift 

fip „ „ 

Aw S at = Wsat - wl = -— (cos fi - A sin /Li) . (4) 

2u>l 


Here A (p) is specific to the type of the topological ob¬ 
ject. For an infinite planar d soliton one has A = 1, cor¬ 
responding to the zero mode on the soliton [21]. In this 


case one finds Aw sa t(/J = 7 t/2 ) = — Au; ma i n (M = 0). For 
a d soliton, separating two HQVs, the signal is modified 
due to finite-size effects and also due to the Aharonov- 
Bohrn effect for spin waves scattered on HQV cores (see 
SM). In this respect, the HQV can be seen as an analogue 
of the Alice string in particle physics [22]. A particle that 
moves around an Alice string flips its charge [23], while 
a quasiparticle going around a HQV flips its spin quan¬ 
tum number. As a result, we expect A < 1. The spin 
polarisation of the HQV core predicted in Ref. 24 does 
not affect the signal, since for two HQVs bound by the 
soliton, the winding of d (and thus spin polarisation) is 
in the opposite directions. 

In a sample, rotating with angular velocity Q, the 
lowest-energy state is achieved when superfluid mimics 
solid-body rotation with an array of rectilinear vortex 
lines, oriented along the rotation direction, with density 
n v = 2Q/(izk). In the polar phase the spin-orbit in¬ 
teraction favours HQVs [25], and they are energetically 
preferable below T c in the axial field (see SM). How¬ 
ever, when solitons are formed between HQV pairs in the 
tilted magnetic field, excess of the spin-orbit interaction 
energy within solitons makes SQVs preferable. Presum¬ 
ably SCVs can be created during cooldown when strong 
time-dependent magnetic field is applied to generate a 
random distribution of vector d [26]. Such experiments 
are beyond the scope of this work. 

We approach the lowest-energy state by slowly cooling 
the sample in rotation from above T c to the polar phase. 
The axis of rotation is parallel to the nafen strands (see 
Fig. la). Results presented in this work were obtained 
using 94% open nafen (density 243 mg/cm 3 ), at 7.1 bar 
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FIG. 2. NMR spectra in the polar phase with HQVs. (a) Nor¬ 
malised spectra measured in transverse field g = 7t/ 2 (blue 
thick line) shows the HQV satellite at the negative frequency 
shift AcOsat and the main line at zero frequency shift. In 
the axial field g = 0 (red thin line) only the main line at 
positive shift Aa; ma i n is seen. This spectrum is not sensi¬ 
tive to presence of vortices, (b) Temperature dependencies of 
the satellite position in the transverse field |Atu sa t(/i = 7r/2)| 
(blue circles) and the main line position in the axial field 
Aoj ma in (g = 0) (red squares) closely match as expected for 
HQVs. The example error bars show full width at half max¬ 
imum of non-Lorentzian main line as an estimate of possible 
systematic error. 


pressure and magnetic field 12 mT. The corresponding 
NMR frequency is ojl/ 27t = 374 kHz. Temperature is 
measured from the NMR spectrum of bulk B phase, using 
the known B-phase Leggett frequency [27, 28]. As a sec¬ 
ondary thermometer we use a quartz tuning fork [29, 30] 
located in bulk 3 He. Details of the experimental setup 
and methods are in the supplementary material [19]. 

In the experiment when the magnetic field is oriented 
transverse to the nafen strands during cooldown, we ob¬ 
serve only the bulk NMR line of the polar phase, which 
is consistent with the creation of SQVs. SQVs do not 
produce a potential well for spin waves and thus yield 
no satellite. When we apply no magnetic field or, alter¬ 
natively, orient the field along the nafen strands during 
cooldown, and then turn the field to the transverse di¬ 
rection below T c , a satellite in the NMR spectrum is ob¬ 
served, Fig. 2a. We interpret this satellite as produced 
by d solitons connecting pairs of HQVs based on the fol¬ 
lowing arguments. 

First, the frequency of the satellite line conforms to 
expectations for a d soliton. In Fig. 2 the frequency 
shift —Aw sat (p, = 7r /2) is compared to Aw ma i n (M = 0) 
as a function of temperature. They are seen to closely 
match each other and we find A (g = 7t/ 2) = 0.93 ± 0.07. 
Deviation from ideal soliton value A = 1 is expected to 
increase with decreasing field angle g due to growing soli¬ 


ton width. This is qualitatively seen in Fig. 3. 

In the analysis we assume that the distance between 
vortex pairs, defining the soliton length, is independent 
of the angle g, i.e., that the tension of the soliton is not 
able to overcome the pinning of HQV hard cores of the 
size of coherence length £ ~ 40 nrn on the nafen strands 
of ~ 10 nrn diameter. The strong pinning is confirmed by 
the observation that after stopping the rotation the satel¬ 
lite in the NMR spectrum remains unchanged for days, 
while the Magnus force, pulling vortices towards sample 
boundary for annihilation, exceeds the soliton tension by 
a factor of 10 3 . 

The second observation in support of HQV is the de¬ 
pendence of the relative satellite peak intensity J sa t on 
the angular velocity fI of rotation, which is applied dur¬ 
ing the cooldown. The measured satellite intensity at 
g = 7t/ 2 is plotted in Fig. 4a. For solitons of the £^> 
width between pairs of vortices, the expected signal is 

_i /q 

/ sa t = (n v /2) g s L^o [9]. Here L = bn v ' is the aver¬ 
age soliton length, where g a ~ 1 is a numerical factor, 
which depends on the distribution of the trapped spin 
wave within the soliton. Numerical factor 6 ~ 1 depends 
on the vortex lattice. For very low vortex density and 
long solitons L —> oo one has g s —> 2. As a result, one 
expects / sat oc Q 1 / 2 , as indeed seen in Fig. 4a. Compari¬ 
son of the measurements with the theoretical prediction 
for J sat [19] shows good quantitative agreement, consid¬ 
ering absence of fitting parameters and simplicity of the 
model, which assumes uniform polar phase and uniform 
distribution of HQVs. 

A remarkable feature seen in Fig. 4a is that the satel¬ 
lite appears in the zero-field cooldowns in the absence of 
rotation as well. We attribute this phenomenon to the 
Kibble-Zurek mechanism (KZ) of vortex (defect) forma¬ 
tion [32, 33] during the crossing of the 2nd order phase 
transition to the polar phase. The KZ mechanism is ex¬ 
pected to create various order-parameter defects includ¬ 
ing vortices of all possible types. In the earlier observa¬ 
tions of the KZ mechanism in superfluid 3 He-B initially 
formed vortices rapidly decay [34-37]. In our case the 
initially formed HQVs freeze due to the strong pinning. 
The scale for the inter-vortex distance is set by the KZ 
length l KZ = £o(tq/t 0 ) 1/4 . Here Tq 1 = - d(T / t Tc) is 


the cooldown rate at T c , £o = £(T = 0) and the order- 
parameter relaxation time To ~ 1 ns. For HQVs the 
inter-vortex distance sets the length of the interconnect¬ 
ing solitons and thus the amplitude of the satellite signal. 
The resulting dependence / sa t oc Ikz oc Tq 1 ^ 4 is indeed 
observed in the experiment, Fig. 4b. The magnitude of 
the signal corresponds to the averaged soliton length of 
1.4 Ikz, as has been estimated as the initial inter-vortex 
distance also in the B phase of 3 He [37, 38]. The shift of 
experimental data in Fig. 4a above the theoretical expec¬ 
tation indicates that the KZ mechanism is important also 
in cooldowns with applied rotation. Detailed analysis of 
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FIG. 3. Frequency shift of the HQV satellite in the NMR 
spectrum, (a) Measured values of the dimensionless frequency 
shift A as a function of the field tilt angle fi (symbols), and 
numerical calculations for the uniform polar phase (solid line) 
using theoretical value of Qj = 17/xm [19]. Leggett frequency 
flp is determined from a separate measurement at fj, = 0. 
As expected, the deviation from the infinitely-long d soliton 
value A = 1 increases towards small fj,. The disagreement 
between the experiment and calculations probably originates 
from disorder in the nafen strand orientation [31], which leads 
to fluctuations of spin-orbit interaction energy within the soli- 
tons. (b) Values of A are found from positions of the HQV 
satellite Aw sa t (blue circles) and of the main line Aa) ma i n (red 
squares). The red and blue solid lines show results of Eqs. (3) 
and (4), respectively, for A = 1. (c) Main line and satellite 
line positions are determined using Lorentzian fits (thick black 
line) of the dispersion signal (red line). For clarity, the lines 
have been shifted vertically. Values of p, are marked above 
the lines. Similar fits are used in extracting the satellite line 
intensity in Fig. 4. In panels (a) and (b) the bars show full 
width at half maximum of the spectral lines. 


the creation of various defects by the KZ mechanism in 
the presence of bias (rotation) and pinning remains a task 
for the future. 

In conclusion, we cool 3 He within nematically ordered 
aerogel-like material called nafen down to the super¬ 
fluid polar phase in rotation up to 2.75 rad/s. When 
the cooldown proceeds without magnetic field or with 
the field oriented parallel to the nafen strands and ro¬ 
tation axis, and afterwards the field is switched on in 
the transverse direction, we observe a satellite peak in 
the continuous-wave NMR spectrum. Dependence of the 
satellite intensity and frequency on the rotation veloc¬ 
ity, temperature, and field orientation identifies it as a 
signal from solitons between pairs of half-quantum vor¬ 
tices. Additionally, the frequency of the satellite reflects 
the Aharonov-Bohm effect for the spin waves scattered by 
the HQV core, which is the analogue of the transfer of the 
Cheshire charge in Alice electrodynamics [22]. The hard¬ 
core HQVs are strongly pinned by the nafen strands, 




FIG. 4. Intensity J sa t of the HQV satellite in the NMR spec¬ 
trum. (a) After cooldown in zero field (red circles) J sa t follows 
expected dependence for HQVs / sa t oc fW 2 . The solid line is 
theoretical prediction for the equilibrium array of HQVs ig¬ 
noring Kibble-Zurek vortices and with no fitted parameters. 
For cooldowns in the transverse field (cyan squares) forma¬ 
tion of HQVs is suppressed. Residual intensity is attributed 
to the possible formation of the soliton sheets terminating at 
the sample walls. The cooldown rate here is Tq 1 = 5T0 -6 s -1 . 
(b) For cooldowns at Q. = 0 rad/s the satellite intensity (sym¬ 
bols) follows dependence 7 sa t oc Tq 1,74 (solid line) expected for 
vortices created by the Kibble-Zurek mechanism, where the 
initial vortex density agrees with earlier measurements in 3 He- 
B [38]. In all cases spectra are measured at /.i = pi/ 2, and 
7 sa t is determined as the area of the satellite normalised to 
the total area of the spectrum. 

which allows us to observe vortices created when passing 
through the transition temperature owing to the Kibble- 
Zurek mechanism even without rotation. If cooldown 
proceeds with the field applied in the transverse direc¬ 
tion favouring stabilisation of single-quantum vortices, 
no satellite is observed as expected. 

In less dense nafen (98% open, 90 mg/cm 3 ) the po¬ 
lar phase acquires axial distortion at low temperatures, 
i.e. it is transformed into the polar distorted A-phase 
[18]. The HQVs there contain unpaired Majorana modes 
in their cores similar to those discussed in HQVs in 
Px + Wy superconductors or in “Kitaev chains” [3-5]. 
Since the HQVs are trapped by the nafen columnar de¬ 
fects, the gap separating the Majorana mode from the 
other Caroli-de Gennes-Matricon excitations localised in 
the vortex core increases from the usual minigap value 
~ A 2 /Ep ~ 10~ 3 A (here Ep is the Fermi energy) in 
an unpinned vortex to a significant fraction of the gap 
amplitude A for a vortex trapped on a defect [39, 40]. 
Our preliminary results of HQV measurements in the 
98% open nafen show that HQVs survive the transi¬ 
tion from the polar phase to the polar-distorted A-phase. 
This opens a pathway for probing the core-bound un¬ 
paired Majorana modes in the experimentally accessible 
temperature range in the polar-distorted A-phase, and 
also in the new topological phases obtained using prop¬ 
erly engineered nanostructured confinement [19], while 
the Kibble-Zurek mechanism makes the HQVs accessible 
with less specialised, stationary equipment. 
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SUPPLEMENTARY MATERIAL 


EXPERIMENTAL SETUP AND METHODS 

The 3 He sample is confined in a 4 x 4 x 4 mm 3 cubic 
container (see Figure la) made from Stycast 1266 epoxy. 
The container volume is filled with aerogel-like material 
called nafen of 243 mg/cm 3 density. The nafen consists 
of strands made from AI2O3, oriented in the same direc¬ 
tion. The average diameter of the strands is measured 
with X-ray scattering as d\ = 9 nrn and variations of the 
diameter are about ±2nm [31]. The average distance 
between strands, determined from the nafen density, is 
c?2 ~ 35 nrn. The nafen was produced by AFN technology 
Ltd in Estonia. The direction of the strands is vertical 
and further denoted as z. The sample is cooled by the 
ROTA nuclear demagnetisation refrigerator [41], which 
can be rotated around its vertical axis. 

Results presented in this work were obtained at 7.1 bar 
pressure and magnetic field 12 mT. The corresponding 
NMR frequency is wl/27t = 374 kHz. The container walls 
and nafen strands have been preplated with a few atomic 
layers of 4 He in order to suppress the paramagnetic signal 
from solid 3 He otherwise formed on all surfaces. At this 
pressure we have measured the transition temperature to 
polar phase T c = (0.9 ± 0.03)T bulk = 1.55 mK, which is 
slightly lower than that reported in Ref. 18 for the same 
nominal nafen density. The difference is probably mostly 
due to variations of the actual density of the nafen. 

The sample is probed using continuous-wave nuclear 
magnetic resonance spectroscopy (cw-NMR). The excita¬ 
tion and signal pick-up is performed with the coils made 
from copper wire which are part of a tuned tank cir¬ 
cuit with a Q value of 150. To improve signal to noise 
ratio a cold preamplifier thermalised to helium bath is 
used. The NMR excitation is applied at the fixed reso¬ 
nance frequency of the tank circuit and the magnitude 
of the magnetic field is swept to record NMR spectra. 
The spectra are averaged over several tens of sweeps in 
order to reduce the noise, except for temperature sweeps 
where averaging is not used. The satellite line position 
and amplitude, and the main line position are found by 
fitting Lorentzian lines to the dispersion signal (see Fig. 
3c), which contains less noise than the absorption. 

Temperature is measured from the NMR spectrum of 
bulk B phase, using known B-phase Leggett frequency. 
The bulk response is picked up by the NMR coils from 
the gaps around the nafen sample and from the filling 
line. It is clearly visible at temperatures above 0.85T C . 
At lower temperatures the bulk line becomes so wide 
that the signal disappears in the noise, and we use as 
a thermometer the position of the NMR line in the po¬ 
lar phase in the axial field. We convert the position to 
the temperature using the model suggested in Ref. 18, 
Eq. (4), where parameters K = 1.2 and T c = 0.9 T c bulk 


are found by comparison to the bulk B-phase tempera¬ 
ture at T > 0.85 T c . As a secondary thermometer we 
use a quartz tuning fork [29, 30] located in bulk 3 He in 
a volume connected by a sintered heat exchanger to the 
nuclear demagnetisation refrigerator. At higher temper¬ 
atures the fork is calibrated from NMR. At the lowest 
temperatures one can also use the onset of the ballis¬ 
tic regime as an approximate reference point [42]. The 
fork we use has the cross-section of a tine 0.4 x 0.3 mm 2 
and the onset takes place at approximately 0.28T C [43]. 
This yields for the lowest temperature reached in the ex¬ 
periments (0.19 ± 0.05) T c , while the lowest temperature 
according to the main NMR line position is (0.2±0.04) T c . 

The magnetic field can be freely rotated in the plane 
perpendicular to the NMR coil axes. This includes field 
orientations parallel to the strands and the rotation axis 
(// = 0), and transverse to them (/j, = 7 t/ 2). The direction 
/i = 7 t/2 corresponds to the maximum frequency shift 
of the satellite Aw sat , see Eq. (4). From this condition 
we have found that the container axis z is tilted by 1° 
from the axis of the magnet assembly. Comparison of 
the NMR spectra at the orientations /1 and fi + tt allows 
us to find the magnitude of the magnetic field trapped 
in the superconducting parts of the setup. Those have 
turned out to be 40 ^T along z, and about 2^T in the 
transverse direction. The container tilt and the trapped 
field are taken into account in the conversion of the NMR 
spectra to the frequency domain. 


ENERGY OF QUANTIZED VORTICES IN THE 
POLAR PHASE 

In a rotating sample the lowest-energy state is achieved 
when superfluid mimics solid-body rotation with an ar¬ 
ray of rectilinear vortex lines. Due to quantized nature of 
vorticity, however, this simulation of solid-body rotation 
cannot be complete. The main contribution to the excess 
energy above the energy of the rotating solid body comes 
from the superflow v s around the vortex core which di¬ 
verges as v s = vk/ (27rr), where r is the distance from the 
vortex axis. If there are spin currents around the core, 
like in the case of HQV, then they also contribute to the 
excess energy. 

The divergence of the superflow is cut at the size of the 
vortex core. In the polar phase of superfluid 3 He both 
SQV (with v = 1) and HQV (with v = 1/2) have hard 
core of the temperature dependent coherence length size 
£(T). Thus their energies per circulation quantum in the 
logarithmic approximation become [44] 

-^SQV = ~^Ps ln j ( S1 ) 
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FIG. SI. A pair of HQVs (green pillars) in the polar phase of superfluid 3 He in transverse magnetic field (blue arrows), 
connected by a soliton in the d field (red arrows). The 7r winding of d around a HQV core is in the opposite direction for the 
two vortices. The soliton width (highlighted by the dash line) is determined by the spin-orbit interaction and is of the order of 
~ 10 p,m, while the intervortex separation ~ 100 pm is set by the rotation velocity applied during the transition through T c . 
Solid lines indicate branch cuts where d —>• — d and the superfluid phase jumps by n, but the order parameter is continuous. 


and 


Ehqy/^ = y + p^)ln|, (S2) 

where p^~ is the transverse component of the superfluid 
density and is the transverse spin rigidity (following 
the notation of Ref. 44). 

The Fermi-liquid effects in superfluid 3 He make p^ < 
Ps, see p. 198 in Ref. 21, with a typical value p^ « 0.5p^~ 
at T = 0.5 T c . Thus in the absence of d solitons, an array 
of HQVs is more energetically favourable than array of 
SQVs, 2Ehqv < Esqv . The ratio p^,/p^ grows with 
increasing temperature and the situation becomes less 
clear when T —> T c and p^, —> p£-. Formally the energies 
of HQV and SQV arrays in the logarithmic approxima¬ 
tion become equal, but the approximation itself ceases 
to be valid, since £ diverges close to T c . To determine 
accurately the energy of vortices one then has to resort 
to numerical calculations of the order-parameter distri¬ 
bution in the core. 

Moreover, since the cooldown in the experiment pro¬ 
ceeds with a finite rate, the vortex formation is not com¬ 
pletely equilibrium process. In this case, to calculate 
the proportion of vortices of each type to be formed, 
one should also take into account the probability to 
form a specific order-parameter configuration in a non¬ 
equilibrium transition through T c . These probabilities 
in general are not directly related to the equilibrium en¬ 
ergies, as was discussed in the context of the transition 


from the normal to the A or B phase of superfluid 3 He 
[45] and can be only estimated in numerical simulations, 
which remain a task for a future research. The clear 
experimental result, though, is that most of the vortices 
formed in the transition to polar phase in rotation in zero 
magnetic field are half-quantum. 

When solitons are formed between HQV pairs in the 
tilted magnetic field, the excess of the spin-orbit inter¬ 
action energy within solitons makes configuration with 
HQVs to be energetically unfavourable at all tempera¬ 
tures. This conclusion is also confirmed by the experi¬ 
ment. 


NMR IN POLAR PHASE OF 3 HE 


The NMR properties of the polar phase, determined by 
the dynamics of the vector d, are similar to that of 3 He- 
A [25]. There are three components of the Hamiltonian 
which are important for spin dynamics: magnetic energy, 
spin-orbit interaction energy, and gradient energy: 

H = F m + F S o + F v , (S3) 


F m = -(S • 7 H) + y XabSaSt, 


Fso = 9d 


A*jA kk + A* k A k j — -A* k Aj k 


F-v = 


Ki{V jA* ak ){S7 jA ak ) 


(54) 

(55) 

(56) 





HQVS AND NMR 


+ KiiVMiVkAaj) + K^VjAljXVkAak )j, 

where S is spin. Susceptibility y is anisotropic, the axis of 
anisotropy is d, and therefore minimum of the magnetic 
energy corresponds to S _L d. 

Let us for simplicity consider the case uniform in the 
z direction. We can parametrise the magnetic and order 
parameter fields: 


m = z, (S7) 

H = h H, h = (y sin/x + zcos/x), 

d = (i cos a + j sin a) sin f3 + h cos /3. 

Here /? is angle between d and the held, and a is the 
azimuthal angle of d in the plane perpendicular to the 
magnetic held. As before, we choose a = 0 along 
the line perpendicular to both H and m which corre¬ 
sponds to the minimum of energy. That is, i = x and 
j = y cos fx — z sin fx. Omitting constant terms, the ener¬ 
gies now become 


Fm = (X-L - X\\)H 2 cos 2 f3, (S8) 

Fso = 2g£>A 2 (sin a sin (3 sin y + cos / 3 cos y) , (S9) 
F v = A 2 K jk [sin 2 /3(Vja)(y k a) (S10) 

+ (Vj/3)(V k p)\, 


where Kj k = K\8j k - 

There are two scales introduced by these energies. The 
ratio of the magnetic and gradient energies gives the mag¬ 
netic healing length £# and the ratio of the spin-orbit and 
gradient energies gives the dipolar length £/> 


& 


2 AT A 2 

H 2 (x± -X||)’ 


AT 
2 9d 


(Sll) 


We use bulk values from Ref. 46 to estimate AT 
scales according to AT oc 1/T C wrt. its bulk value (see 
Ref. 47), and ignoring the additional scattering due to 
the nafen we get £d = 17 /im. Note that our dehnition 
of follows the notation of Ref. [44], which is smaller 
by y/2 than that in e.g. Ref. [46]. 

In the high-field limit £,d ^ £h, that is, magnetic en¬ 
ergy is in the minimum everywhere excluding small re¬ 
gions of the £h size (for example, cores of spin vortices). 
The small volume of this regions makes them invisible in 
NMR experiments. In the rest of the volume f3 = 7r/2, 
and only variations of a are important. Minimising the 
total energy wrt. a yields Eq. (2). In a tilted magnetic 
held there are two possible uniform equilibrium textures 
with a = 0 and a = n, that is, vector d is perpendicular 
to both H and 1 and can point in two possible direc¬ 
tions. The two states may coexist, and in between there 
is a quasi-stable d soliton. The thickness of such soli- 
ton can be solved analytically and is of the order of £r>. 
The soliton may terminate e.g. at a container wall, at a 
half-quantum vortex, or at a spin vortex. 


The NMR response of bulk polar phase can be derived 
by considering tiny oscillations of the spin part of the 
order parameter held. The position of the satellite peak, 
derived in a similar fashion, is given by the eigenvalue 
A of the wave equation describing the spin wave modes 
localised in the texture of the vector d (see Fig. SI). Note 
that we dehne A so that it is positive for practical reasons. 
Assuming vortices are oriented along z, it acquires the 
following form [44]: 


U( r) 




4/ ,(S12) 


—Ci)(V Q ') 2 — sin 2 asin 2 y , A(r) = Va ,(S13) 

H 2 

w sa t -co L = —^(-Asin 2 /x + cos 2 /x) (S14) 
2 co L 


For a pure inhnite soliton in the (y, z)-plane one has 
a(x = +oo) — a(x = —oo) = ir. In this arrangement 
the vector potential A(x) = Va(x) can be gauged away 
by the transformation 4/ = Then one obtains 

A = 1, which rehects the mathematical connection of 
the localised mode to another collective bosonic mode 
trapped in the soft core of the soliton — the zero mode 
corresponding to the translation of the soliton with the 
eigenfunction 4/(a;) = da/dx [21]. 

In the ^D-vicinity of a half-quantum vortex a ss <p/2+ 
const. (<p is the azimuthal coordinate). Therefore, po¬ 
tential U(r < £d) oc l/(4r 2 ), where r denotes the dis¬ 
tance from the HQV core. This would lead to a ground 
state much more shifted from the bulk line than the one 
observed, and the actual position would depend on the 
precise way U is cutoff at r ~ £ inside the vortex core. 
This large shift of the order of is, however, com¬ 

pensated by the so-called Aharonov-Bohm effect, which 
is manifested by the ir phase jump of the wave func¬ 
tion around the vortex core. The analogy to the original 
Aharonov-Bohm effect can be seen by introducing an ef¬ 
fective magnetic field concentrated in the vortex cores 
[44]: 


B c ff = Vx A = 27wzc>(r). (S15) 

Consequential additional scattering for spin waves leads 
to a reduced shift of the satellite line in the NMR spec¬ 
trum. 

In general, the satellite peak position depends on the 
aspect ratio of the soliton L/D, where D is the thickness 
of the soliton. In the tilted magnetic field the thickness 
of the soliton scales as D ~ £d/ sin /x and thus the finite 
size effect becomes important at smaller /x. We calcu¬ 
lated A numerically as a function of L sin /x/^/j in a 2D 
finite element model, i.e., assuming infinite straight vor¬ 
tices and uniform polar phase. Five different cases were 
studied (Fig. S2). Periodic structures of infinite soliton 
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a 


d 



L sin |_i/^ D 


L sin |_i/^ D 


FIG. S2. Calculated values of g s (bottom left) and A (bottom 
right) as a function of L sin g/^D ■ We consider five cases: (a) a 
single soliton of a finite length L, periodic structures of infinite 
solitons layers with the period L and same (a) or alternating 
soliton orientations (b), and square lattices of solitons and 
half-quantum vortices with periodic boundaries such that the 
soliton orientations are alternating (d) and not alternating 
(e). 


layers yield A = 1, as expected for any unpinned soliton 
structure. Periodic square-lattice structures of solitons 
terminating to half-quantum vortices, and a single soli¬ 
ton between two HQ Vs all yield roughly the same A (//). 
Details of the calculation will be published separately 

[48]. 

The relative intensity of the satellite peak is obtained 
using the eigenfunction of the bound state in the soliton 
(see Eqs. (13) and (39) in Ref. 9). We get 


Fsat — - tl 


1 |/ ^>d 2 r | 


f2 

■=D 


f |\I/| 2 d 2 r 


2 sin 2 /i 


n v I 


-M 


(S16) 


Here I M = g s L sin/.</££> where the numerical factor g s 
depends on, e.g., the shape of the wave function. The 
equilibrium density of vortices is 


n v 


2TO3H 

virh 


(S17) 


The distance L between the HQVs is given by L 2 = 
b 2 /n v , where 6=1 for a square lattice used in the simula¬ 
tions. In a transverse field the relative satellite intensity 
from the HQVs becomes 


I s at = ( S18 ) 

For = 17 /im our numerical simulation yields g s ss 1.7 
independent of the configuration. 


TOPOLOGY AND SUPERFLUID PHASES OF 
3 HE 


Many exotic topological objects have been observed 
experimentally in the superfluid phases of 3 He, among 
them quantized vortices with spontaneous asymmetry 
in the vortex core, skyrmions with different topological 
numbers, soliton terminating on the combined spin-mass 
vortices, and vortex sheets of different geometry. Some 
exotic objects are still to be detected, such as a nexus - 
a Dirac monopole terminating a vortex line, or a boojum 
- topological defect living on the surface of the container 
or at the interface between different superfluid phases [1]. 

The known phases of liquid 3 He belong to 4 different 
topological classes: 

(i) The normal liquid 3 He belongs to the class of sys¬ 
tems with topologically protected Fermi surfaces. The 
Fermi surface is described by the first odd Chern num¬ 
ber in terms of the Green’s function [1]. 

(ii) Superfluid 3 He-A and 3 He-Ai are chiral superfluids 
with topologically protected Majorana-Weyl fermions in 
the bulk. In the relativistic quantum field theories, Weyl 
fermions give rise to the effect of chiral (axial) anomaly. 
The direct analogue of this effect has been experimen¬ 
tally demonstrated in 3 He-A [49]. It is the first condensed 
matter system where the chiral anomaly effect has been 
observed. The singly quantized vortices in superfluids 
with Weyl points contain dispersionless band (flat band) 
of Andreev-Majorana fermions in their cores [50]. Dis¬ 
coveries of Weyl fermions have been recently reported in 
the topological semiconductors [51, 52], 

(iii) 3 He-B is the purest example of a fully gapped 
Dill class superfluid with topologically protected gapless 
Dirac fermions on the surface [53]. Some evidence of gap¬ 
less fermions on the surface of 3 He-B has been reported 
[54, 55], 

(iv) Polar phase, the most recently discovered phase, 
belongs to the class of fermionic materials with topolog¬ 
ically protected lines of nodes, and thus two-dimensional 
flat bands of Andreev-Majorana fermions cover its sur¬ 
faces [56]. (See also recent reviews on superconductors 
with topologically protected nodes [57].) Recently the 
class of these materials has been joined by the topologi¬ 
cal semimetals [58-60]. The half quantum vortex, which 
we found in the polar phase, can be represented as a vor¬ 
tex in a single spin component. The spectrum of bound 
states experiences condensation of zero modes [61]. 

It is possible that with properly engineered nanostruc- 
tured confinement one may reach also new topological 
phases of liquid 3 He including 

(i) The planar phase, which is non-chiral, has Dirac 
nodes in the bulk, and Fermi arc of Andreev-Majorana 
fermions on the surfaces [62]. 

(ii) The two-dimensional topological states in the 
ultra-thin film, including inhomogeneous phases of su¬ 
perfluid 3 He films [63]. The films with the 3 He-A and 
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the planar phase order parameters belong to the 2D fully 
gapped topological materials, which experience the quan¬ 
tum Hall effect and the spin quantum Hall effect in the 
absence of magnetic field [64]. 

(iii) a-state, which contains 4 left-handed and 4 right- 
handed Weyl points in the vertices of a cube [65]. It 
resembles the high energy physics model system with 8 
left-handed and 8 right-handed Weyl fermions in the ver¬ 
tices of a four-dimensional cube [66]. 
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